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DESTRUCTION OF STELLAR DISKS BY PHOTOEVAPORATION 
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RESUMEN 

El resumen sera traducido al espaiiol por los editores. Photoevaporation may provide an explanation 
for the short hfetimes of disks around young stars. With the exception of neutral oxygen lines, the observed 
low- velocity forbidden line emission from T Tauri stars can be reproduced by photoevaporating models. The 
natural formation of a gap in the disk at several AU due to photoevaporation and viscous spreading provides 
a possible halting mechanism for migrating planets and an explanation for the abundance of observed planets 
at these radii. 

ABSTRACT 

Photoevaporation may provide an explanation for the short lifetimes of disks around young stars. With the 
exception of neutral oxygen lines, the observed low-velocity forbidden line emission from T Tauri stars can 
be reproduced by photoevaporating models. The natural formation of a gap in the disk at several AU due to 
photoevaporation and viscous spreading provides a possible halting mechanism for migrating planets and an 
explanation for the abundance of observed planets at these radii. 
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1. INTRODUCTION 

Most low-mass stars form surrounded by circum- 
stellar disks. Observations of infrared excess (e.g., 
Haisch, Lada, & Lada 2001) suggest disk lifetimes 
Tdisk ~ 6 X 10^ yr, requiring an efficient mechanism 
for disk dispersal. 

HoUcnbach, Yorke, & Johnstone (2000) consid- 
ered a variety of disk dispersal mechanisms and con- 
cluded that viscous accretion of the disk onto the 
central star (e.g. Hartmann et al. 1998) together 
with photoevaporation of the disk at moderate radii 
(Shu, Johnstone, & HoUenbach 1993; HoUenbach et 
al. 1994; Johnstone, HoUenbach, & Bally 1998) must 
act together efficiently to remove the entire disk. 
Numerical calculations by Clarke, Gendrin, & So- 
tomayor (2001) and Matsuyama, Johnstone, & Hart- 
mann (2003a) have shown that these processes may 
remove the disk in 10^"^ yr. 

The photoevaporation model fits observational 
data well in the case of external heating via nearby 
massive stars (Bally et al. 1998; Johnstone, HoUen- 
bach, & Bally 1998; Storzer & HoUenbach 1998). 
With the exception of a few cases (e.g., MWC349A), 
evidence for disk photoevaporation due to the cen- 
tral star (Shu et al. 1993; HoUenbach et al. 1994) is 
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largely circumstantial. Observations of blue-shifted, 
low-velocity emission from forbidden lines of oxy- 
gen, nitrogen, and sulfur in the spectra of many T 
Tauri stars (Hartigan, Edwards, & Ghandour 1995) 
provide useful diagnostics. Font et al. (2004) cal- 
culated the flow properties of the photoevaporative 
disk wind and found them to compare reasonably 
with the strengths and profiles of the nitrogen and 
sulfur lines. The oxygen lines, however, are under- 
abundant in the model. 

Along with disk dispersal, photoevaporation and 
viscous accretion produce structure in the disk such 
as gaps and rings. Matsuyama, Johnstone, & Mur- 
ray (2003b) showed that the formation of gaps within 
the gaseous disk during the dispersal era places con- 
straints on the migration of planetary orbits. 

The following sections review the key results of 
Matsuyama et al. (2003a,b) and Font et al. (2004). 

2. PHOTOEVAPORATING DISKS 

The evolution of a stellar disk under the influence 
of viscous accretion and photoevaporation from the 
central source or external stars has been studied by 
Clarke et al. (2001) and Matsuyama et al. (2003a). 
The key parameter in these studies is the source and 
strength of the ionizing radiation Clarke et al. 
assumed that the ionizing flux was due to super-solar 
activity on the surface of the star and constant with 
time. Matsuyama et al., however, treated the pho- 
toionizing flux from the central source as arising from 
both the quiescent star and accretion shocks at the 
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base of stellar magnetospheric columns. In the lat- 
ter study, therefore, the ionizing flux was calculated 
self-consistently from the accretion mass-loss rate. 

Alexander, Clarke, & Pringle (2004) raised objec- 
tions to the use of accretion energy to produce ultra- 
violet photons, showing that such photons would be 
absorbed in the dense accretion flow and may not es- 
cape the shock front. Furthermore, both stellar and 
accretion shock ultraviolet photons may be trapped 
by the strong protostellar jet which launches from 
the inner disk. Shang et al. (2002) were unable to 
ionize completely the jet using only photons from 
the central star due to the high optical depth of the 
jet itself. It may be that the source of ultraviolet 
radiation ionizing the disk is the disk-jet boundary, 
in which case the photon flux still scales as the ac- 
cretion energy. Other possibilities include enhanced 
stellar surface activity (but only in stars with weak 
or absent jets) or nearby massive stars. 

The disk cannot be entirely removed using only 
viscous accretion and photoionization from the disk- 
star accretion shock because the photoionizing flux 
decreases too quickly with time (Matsuyama et al. 
2003a). When 6-13.6 eV far-ultraviolet (FUV) pho- 
tons from relatively nearby massive stars can be in- 
cluded, however, the disk is removed in lO^^^yr. In 
addition, when ionizing photons from close massive 
stars can be included the disk is removed in 10^~^ yr. 

An intriguing consequence of photocvaporation 
by the central star is the formation of a gap in 
the disk at late stages of the disk evolution (Figure 
The gap forms near rg ~ GMi,/a?, where the 
sound speed a of the 10^ K photoevaporated disk gas 
matches the escape velocity from the central star. In- 
terior to this radius, material from the disk is bound 
to the system and forms a hot disk corona. Exterior 
to this radius, the ionizing radiation is less intense 
due to divergence. The gap forms when the surface 
density at drops to S(rg) ^ 2gcm~^, for typical 
ionization and viscous parameters. 

When viscous accretion and photocvaporation by 
both the central star and nearby massive stars are 
considered simultaneously, the disk shrinks and is 
truncated at an outer radius, set by the temperature 
in the FUV photoevaporated disk gas (T 10^ K) 
or the temperature in the ionizing photon photoe- 
vaporated flow [T ^ 10"^ K). At typical distances 
{d ^ 0.03 pc) from massive stars the disk is evap- 
orated by both external FUV photons and internal 
ionizing photons. This produces both an outer limit 
for the disk size and an inner disk gap during the 
last stages of disk evolution (Figure [21) . 
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Fig. 1. Snapshots of the disk surface density for a repre- 
sentative model with both viscous evolution and photoc- 
vaporation from the central star. Rgii refers the gravita- 
tional radius for 10'' K gas. See Matsuyama et al. 2003a 
for details. 
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Fig. 2. Snapshots of the disk surface density with vis- 
cous evolution, internal photocvaporation, and evapora- 
tion due to a nearby massive star. Rgi refers the gravita- 
tional radius for 10^ K gas. See Matsuyama et al. 2003a 
for details. 

3. HYDRODYNAMIC MODELS 

Low-velocity (e.g. ~ lOkms^^) blue-shifted for- 
bidden lines are observed in the spectra of T Tauri 
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Fig. 3. Steady-state flow results for the photoevaporative 
disk wind model. The solid lines are streamlines and 
indicate how the ionized gas travels. The dashed lines 
indicate where vtot/a equals 1, 2, and 3. See Font et al. 
(2004) for details. 

stars (Hartigan et al. 1995). The similarity between 
the blue-shift of the low- velocity component and the 
sound speed in 10'' K gas has led to speculation that 
this feature may be an observational signature of 
photoionized disk winds. 

Using the analytic model of Shu, Johnstone, & 
HoUenbach (1993) and HoUenbach et al. (1994) as 
a basis, Font et al. (2004) examined the character- 
istics of photoevaporative outflows using hydrody- 
namic simulations. Figurc|31 shows the streamlines 
and isotachs of the resultant disk wind. Near the 
disk surface, the radial density gradient forces the 
streamlines to bend outward significantly and the 
flow accelerates. At larger distances from the disk, 
however, the flow becomes radial and approaches the 
classical Parker wind solution (Parker 1963). 

General results from the simulations were found 
to agree well with the analytic predictions, although 
some small differences were present. Most impor- 
tantly, the flow of material from the disk surface 
develops at a somewhat smaller radius than in the 
analytic approximations {rg = G M^/'Ha?) and the 
flow-velocity from the disk surface is only one-third 
the sound speed. 

Assuming the gas to be almost entirely ionized 
and adding contributions from each cell in the hy- 
drodynamic simulation, Font et al. (2004) computed 
the widths, velocities, and luminosities of sulfur and 
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Fig. 4. Comparison of [SII]A6731 line profiles. The pho- 
toevaporative disk wind (PDW) model profile is calcu- 
lated assuming 0, — 10^^ s~^, M* — 1 Mq and a face-on 
disk. See Font et al. (2004) for details. 

nitrogen forbidden lines. The simulated line shapes 
were found to be in relatively good agreement with 
observations (Figure ^ . 

Font et al. (2004) also demonstrated that the 
model line strengths are in agreement with ob- 
served low- velocity forbidden line emission of ionized 
species from T Tauri stars. This is in contrast with 
magnetic wind models (Garcia et al. 2001a, 2001b), 
which systematically under-predict these line lumi- 
nosities. The photoevaporative model, however, can- 
not account for the luminosities of neutral oxygen 
lines in T Tauri stars since almost all oxygen in the 
model is ionized. 

4. MIGRATION OF PLANETARY ORBITS 

The recent discovery of Jupiter-mass planets or- 
biting at 3-5 AU from their stars complements earlier 
detections of massive planets with very small orbits. 
The short period orbits suggest that migration has 
occurred, possibly due to tidal interactions between 
the planets and the gas disk. The newly discovered 
long period planets (and our own solar system) show 
that migration is either absent or rapidly halted in 
at least some systems. 

Matsuyama et al. (2003b) proposed a novel mech- 
anism for halting disk coupled planet migration at 
several AU in a gas disk. As mentioned in ^ pho- 
toevaporation of the disk may produce a gap in the 
disk at rg (e.g. a few AU). Such a gap would pre- 
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Fig. 5. Planet streamlines in a viscously evolving and 
photoevaporating disk. See Matsuyama et al. (2003b) 
for details. 

vent outer planets from migrating inward toward the 
star, resulting in an excess of systems with planets at 
or just outside Tg. Figure [3 plots the streamlines of 
planets moving in tandem with their viscously evolv- 
ing disk. Planets which form in the inner disk mi- 
grate with the viscous disk toward the central star. 
In contrast, planets from the outer disk are halted 
by the gap near r^, fixing their orbits at several AU. 
Almost all planets which tidally lock themselves to 
the disk early will first migrate outward in the disk 
before viscously accreting toward the star. 

5. CONCLUSIONS 

Photoevaporation of disks around young stars 
may be responsible for short observed disk lifetimes, 
especially when powered by FUV or ionizing radi- 
ation from nearby massive stars. Together, pho- 
toevaporation and viscous accretion naturally lead 
to the formation of gaps and ring structures within 
disks, without the need to invoke unseen planets. 

Hydrodynamic models of the photoevaporative 
disk wind reveal that the launching point for the 
flow, Tg, has been overestimated in the analytic mod- 
els by a factor of ~ 3 and that the launch velocity 
has been also somewhat overestimated. While most 
observed low-velocity forbidden line profiles are rea- 
sonably matched by the model, the predicted neutral 
oxygen forbidden line Hux is too low since almost all 
oxygen in the model is ionized. Interestingly, a hot 
10^ K wind which is not photoionized would contain 



mostly neutral oxygen and the line strengths would 
produce an excellent fit to the observations. 

The existence of a gap in the disk at Tg provides a 
halting mechanism for migrating planets. Given the 
predicted location of Vg at several AU from the cen- 
tral star, such gaps might explain the recent abun- 
dance of observed planets with these radii. 
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